Determinants related to Dirichlet characters 
modulo 2, 4 and 8 of binomial coefficients and the 
algebra of recurrence matrices 
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| Abstract: Using recurrence matrices, defined and described with some 

details, we study a few determinants related to evaluations of binomial co- 
efficients on Dirichlet characters modulo 2, 4 and 8. N 

i—l . 

1 Introduction 

The aim of this paper is twofold: It contains computations of a few de- 
terminants related to binomial coefficients. The most interesting example, 
discovered after browsing through [7], is obtained by considering binomial 
coefficients modulo 4. We include two similar examples taken from [2] and 

The second topic discussed in this work are recurrence matrices, see [3] 
for a very condensed outline. They are defined as certain sequences of matri- 
ces involving self-similar structures and they form an algebra. Computations 
in the algebra of recurrence matrices are the main tool for proving the de- 
terminant formulae mentionned above. Recurrence matrices are however of 

00 ' 

■ independent interest since they are closely linked for example to automatic 

C > . sequences, see PQ, to rational formal power series in free noncommuting 

variables and to groups of automata, see [S] for the perhaps most important 
example. The determinant calculations of this paper can thus be considered 
as illustrations of some interesting features displayed by recurrence matrices. 
H ' The sequel of this paper is organised as follows: 

The next section recalls mostly well-known facts concerning binomial 
and q— binomial coefficients and states the main results. 

Section [3] defines the algebra 1Z of recurrence matrices. It describes them 
with more details than necessary for proving the formulae of Section [2j 
Section [J] discusses a few features of the group of invertible elements in 

K. 
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Section [5] proves formulae for the determinant of the reduction modulo 
2 of the symmetric Pascal matrix (already contained in [2]) and of a deter- 
minant related to the 2— valuation of the binomial coefficients (essentially 
contained in [3]). 

Section [6] is devoted to the proof of our main result, a formula for 
det(Z(n)) where Z(n) is the matrix with coefficients X-b(( s ^*)) £ {0, ±1}, < 
s,t < n obtained by considering the "Beeblebrox reduction" (given by the 
Dirichlet character XB(^ m ) = 0, XB(4m ± 1) = ±1 for m E Z) of binomial 
coefficients. The proof uses an LU factorisation of the infinite symmet- 
ric matrix Z = Z{oo) and suggests to consider two (perhaps interesting) 
groups and Tz whose generators display beautiful "self-similar" struc- 
tures. This section ends with a short digression on the "lower triangular 
Beeblebrox matrix" and the associated group. 

Section [7] contains some data concerning the reduction of binomial coef- 
ficients by a Dirichlet character modulo 8 related to the Jacobi symbol. 

Section [8] reproves a known formula for evaluating q— binomial coeffi- 
cients at roots of unity. This formula yields easily formulae for some deter- 
minants associated to the reduction modulo 2 and the Beeblebrox reduction 
of (real and imaginary parts) of q— binomial coefficients evaluated at q = — 1 
and q = i. 



2 Main results 

2.1 Reductions modulo 2 

Let P{n) be the integral symmetric n x n matrix with coefficients P s> t 6 
{0, f }, < s, t < n defined by 

P s ,t^( S ^) (mod 2) 

where ( s ^*) = ^~fr denotes the usual binomial coefficient involved in the 
expansion (x + y) n = £Lo {l)x k y n ~ k . 

The evaluation ( s ^*) (mod 2) can be computed using a Theorem of Lu- 
cas, see [TT], page 52. Given a prime number p, it states that 



Y[ ( 3 ) (mod p) 
j>0 ^ 



where z/j, m £ {0, 1, ... ,p — 1} are the coefficients of the p— ary expansion of 
n = ^jx^Vjp 3 and k = ^2j>o KjP 1 ■ Another formula (due to Kummer) for 
(k) (mod 2) will be presented in Section 12.21 

Let ds(n) = £}=o 2(n)J Vj 6 N denote the digit -sum of a natural integer 
with binary expansion n = J2j>o v fl i \ u o-> v \i ' ' ' £ {0> !}• 
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Theorem 2.1. We have 



det(P(2n)) = (-1)" 

and 

det(P(2n+ 1)) = (-i) n +Mn) _ 

Remark 2.2. T/ie infinite symmetric integral matrix P with coefficients 
Ps,t = ) given by the binomial coefficients is sometimes called the Fer- 
mat matrix. Vandermonde's identity ^2 k=0 = X^=o (k) (t-fc) = (^f) 

shows that det(P(n)) = 1 where P(n) is the symmetric nxn submatrix with 
coefficients ( s ^*) , < s,t <n of P. 

2.2 2— valuations 

Given a prime p, we denote by v p : Q* — > N the p— valuation. Any rational 
number a can thus be written in the form a = p Vp ^^ with n,m £ Z 
coprime to p. Let V(n) be the symmetric nxn matrix with coefficients 
V s>t G {±l,±i} given by 

V Sit = i v < s V)), 0< s ,t<n . 

The p— valuation ^p(( s ^*)) of a binomial coefficient can be computed 
using a Theorem of Kummer stating that equals the number of 

carries occuring during the addition of the p— ary integers s = 2^j>o a jP J 
and t = Ylj>o T jP J ■ More precisely, Kummer shows the identity 

(see Lehrsatz, page 115 of [10]) where x = YljXj^jP 3 e N with £j £ 
{0, 1, ... ,p — 1}. This implies the formula 

Mi S+ j) = v p ((s + t)\) -Vp(sl) -v p (t\) = —^—^(crj+Tj-Uj) 
V s J p j>o 

(see [10], page 116) where aj,Tj,Uj £ {0, 1, ... ,p — 1} are defined by the 
p-ary expansions s = Ylj>o a jP 3 \t = Ej>o T iP j and s + t = J2j>o u jP j ■ 

The next result uses the regular folding sequence / : {1, 2, . . . } — > {±1}- 
It is defined recursively by /(2 n ) = 1 and f(2 n + a) = —f(2 n — a) for 
1 < a < 2™, see for example pQ. 

Theorem 2.3. We have 

2n-l 

det(V(2n)) = (-1)" [] (1 - f(k)i) G Z[i] 
fc=i 
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and 

2n 

det(V(2n + 1)) = JJ(1 - f(k)i) G Z[i] 

fc=i 

(twit/i ds(^2j =Q Vj2 J ) = X^j>o denoting the binary digit-sum). 

Remark 2.4. Xei Z? denote the diagonal matrix with diagonal entries i ds (°\ i ds ^\ i ds ( 2 \ 
The paper J3]/ deals with the Hankel matrix H defined by H s j = i ds ( s+t ) , < 
s, t related by H = DVD to the complex conjugate V of the matrix V in- 
volved in Theorem \2.3[ 

Let us also mention that slight extensions of the computations occuring in 
our proof of Theorem \2.3\ establish the existence of nice continued J— fraction 
expansions for the formal power series (cf. J$) 

oo 
fc=0 

1 / 1 + i 1 — X f-r , 9 & . \ 

i — +— n< i +- 2 ) . 

\ k=0 / 

1 / 1 + i i — 1 1 — x 2 j-j . . 2 fe. \ 

\ k=0 / 

2.3 Beeblebrox reduction 

The idea (and the "Beeblebrox" terminology) of considering the "Beeblebrox 
reduction" of binomial coefficients are due to Granville, see [7J and [8]. 

We define the "Beeblebrox reduction" as the Dirichlet character \B '■ 
Z — ► {0, ±1} given by 

{0 if x = (mod 2) 
1 if x = 1 (mod 4) 
-1 if a = 3 (mod 4) 

or equivalently by xb(2Z) = 0,x_b(4Z ± 1) = ±1. Beeblebrox reduction is 
the unique Dirichlet character modulo 4 not factorising through Z/2Z. 
The following result allows fast computations of \B ( \T\ ) • 

Theorem 2.5. We have 



(_l)n( fc+ l) XB( (n )} 

We denote by Z'(re) (where the letter Z stands for Zaphod Beeblebrox, 
following the amusing terminology of [7J and [8]) the symmetric Beeblebrox 



x*(( 2 2 p) = 

Xs(( 2 fc+i)) = 
X*(GK)) = 
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matrix of size n x n with coefficients Z s j G {—1, 0, 1} for < s, t < n given 
by the Beeblebrox reduction Z s j = Xs(( s+ *)) °f binomial coefficients. 
Define / : N — > ±3 Z by /(0) = 1, /(l) = -1 and recursively by 



f f2 a + 6) - I 3/(6) if 26 < 2 ° 
71 J 1 lf(b) otherwise 



for n = 2 a + b > 2 where < 6 < 2 a . Laurent Bartholdi pointed out that 
the value f(n) of a binary integer n = J2i>o ^ s a ^ so gi ven by 

f( n ^ = I Ui=0, v i+1 =l}-${i | Ui=u i+1 =l} 

= (-l) n 3( 1 -2 Jy i)^+l 
i>0 

Theorem 2.6. M^e /icwe 

n-l 



det(Z(n)) = J] f(k) G ±3 N . 



fc=0 



2.4 The Jacobi-character modulo 8 

Let xj '■ N — ► {0, ±1} denote the Dirichlet character modulo 8 defined by 
Xj(2Z) = and 

f 1 if n = ±l (mod8), 

XJl j \ -1 if n = ±3 (mod 8) . 

Since xj(p) = 2^ 1 )/ 2 (mod 2) for p an odd prime, we call xj the Jacobi- 
character and we consider the matrix J(n) with coefficients J s j = Xj(( s ^*)) 
for < s, t < n. The techniques of this paper can be used to prove the 
following result: 

Theorem 2.7. We have 

n-l 

det(J(n)) = Yl g{k) G ±3 N 
fc=0 

where 

g{n) = (-l) n Y[r ([n/2kl) 

k>0 

with e(k) the 8— periodic function given by the table 



k 


1 2 3 4 5 6 7 


e(k) 


1 -1 -2 3 -1 
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The matrices Z(n) and J(n) share many features. However the consid- 
erable greater complexity of all objects attached to J makes it more difficult 
to pin down interesting algebraic structures associated to J. 

Let me also add that the remaining Dirichlet character modulo 8 (given 
by x(2Z) = and x(ea + 8Z) = e for a G {1,3} and e G {±1}) does not 
seem to give something interesting: the symmetric infinite matrix obtained 
by applying % to binomial coefficients has probably no LU decomposition 
in the algebra of recurrence matrices. 

Similarly, there seem to be no new interesting Dirichlet characters 
(mod 16) (with values in the Gaussian integers {±1,±\/— T}- 

Remark 2.8. The case of the character n i — ► x{ n ) = e n = n^ 1 ^ 2 
(mod p) with e n G {0, ±1} /or p an odd prime gives rise to similar re- 
sults which are somewhat trivial. Indeed, Lucas's theorem implies that the 
corresponding infinite matrix with coefficients x(( l ~^)), — J has a struc- 
ture of an infinite tensor-power ( corresponding to a "recurrence matrix of 
complexity 1"). It is thus not very interesting and easy to handle. 

The same remark holds for the remaining characters (mod p) when 
working in a suitable field (or integral subring) of cyclotomic numbers. 



2.5 q— binomials 

The expansion (x + y) n = Ylk=o © x k y n ~ k involving two non-commuting 
variables x, y related by yx = qxy where q is a central variable defines the 
q— binomials coefficients 

J =713— ^-wZTTT- 



(nj=i(i -<?')) (n-=i fc (i-^)) 



An ordinary binary coefficient ( s ^*) can be identified with the number of 
lattice paths with steps (1,0) and (0,-1), starting at (0, s) and ending at 
(t, 0). Similarly, the coefficient of q c in the q— binomial ( s ^"*) ? counts the 
number of such paths delimiting a polygon of area c in the first quadrant 
{(x,y) £K 2 | x,y > 0}. 

Reflecting all paths contributing to ( s ^*) 9 with respect to the diagonal 
line x = y yields the equality 

s + _ ^s + t 

Rotating all paths contributing to ( s ^ t ) Q by a half-turn centered at \(t, s) 
shows the identity 

s + t\ st fs + t 



Q 

s 
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their first step (horizontal or vertical) shows the recursive formula 



Partitioning all paths contributing to ( s ]~ ) accordingly to the nature of 



or equivalently (*£) q = ? fe ( n fc ) 9 (fe-i) which is the q— version of the 
celebrated recurrence relation (?) = ( n ^. 1 ) + (uZi) f° r ordinary binomial 
coefficients. 

Cutting all lattice paths 7 contributing to ) along the diagonal line 
s = t in two lattice paths shows the formula 

where k E {0, 1, . . . ,min(s,t)}. This identity amounts to the matrix iden- 
tity P q = L q D q L q where P q is the infinite symmetric matrix with coef- 
ficients ( s ^ t ) q , < s,t, where L q is the lower triangular unipotent ma- 
trix with coefficients (*) , < s,t and where D q is diagonal with diag- 
onal coefficients 1, q, g 4 , g 9 , g 16 , g 25 , .... Denoting by P q (n) the submatrix 

( S ^s *) ' ^ — s > ^ < n formed by the first n rows and columns of P q we 

™-i -2 

have the identity det(P g (n)) = q^3= a 3 which specialises to the identity 
det(Pi(n)) = 1 of Remark 12.21 Appendix I of [6] contains many more for- 
mulae for (?J . 

The formula of Lucas 



a \ _ fl a /p\\ f a (modp)\ 
b) = {[b/ P ]) [b (modp)J {m ° dp) 

(where p is a prime number and where a (mod p), 6 (mod p) E {0, 1, . . . ,p — 
1}), see Section [2~T1 or [11J . has the following known analogue for g— binomials 
which reduces their evaluation at roots of 1 of small order to evaluations of 
ordinary binomial coefficients. 

Theorem 2.9. If uj = e 2lnk / n is a primitive n—th root of 1 (ie. (k,n) = 1 
with k E Z and n E N ) then 

/ a\ f[a/n\\ (a (mod n) 
b) ul = \[b/n\) 1 \b (modn) 

for all a, b E N where a (mod n), b (mod n) E {0, 1, . . . , n — 1}. 



Theorem 12.91 can be used to establish formulae for determinants of the 
symmetric matrices obtained by considering the reduction modulo 2, the 
Beeblebrox reduction or the reduction using the Jacobi character modulo 8 
of (the real and imaginary part of) ) , < s, t < n evaluated at q = — 1 
and q = i. 
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3 The algebra of recurrence matrices 



Recurrence matrices, introduced in [3], are a convenient tool for proving 
our main results. Recurrence matrices are closely related to rational formal 
power series in free non-commutative variables and can be considered as 
generalisations of finite state automata or of iterated tensor products. They 
arise also naturally in the context of "automata groups" , a notion generalis- 
ing a famous group of Grigorchuk, see [§] . The following exposition does not 
strive for exhaustivity or for the largest possible generality. Generalisations 
(e.g. by replacing the field of complex numbers by an arbitrary commutative 
field or by considering sequences of square matrices of size k n x k n ,n E N 
for k £ {1, 2, 3, . . . }) are fairly straightforward and contained in [3] or with 
more details in [3]. 

The papers [5] and [12] deal with interesting subalgebras, called self- 
similar algebras, formed by recurrence matrices. 



3.1 Recurrence matrices 

Consider the vector space 



A= [J M 2«X2"(C) 



n=0 

whose elements are sequences A = (A[0], A[l], A[2], . . . ) with A[n] £ M2«x2 n (C) 
denoting a complex square matrix of size 2 n x 2 n . The obvious product 

AB = (A[0)B[0],A[1)B[1],A[2)B[2],...) 

turns A into an associative algebra. Denoting by 

p{0, 0) A p(0, 1)A, p(l, 0)A, p(l, I) A e A 

the four "corners" of 

(p(0,0)A)[0] (p(0A)A)[0] \ f (p(0,0)A)[l] (p(0,l)A)[l] ' 



A A[0] >\ [p(l,0)A)[0] (p(l,l)A)[0] J>\ (p(l,0)A)[l] (p(l,l)A)[l] 

obtained (after deletion of the lxl matrix A[0]) by considering for all n > 1 
the 2 n_1 x 2 n ~ 1 submatrix defined by the first or last 2 n_1 rows and by 
the first or last 2 n_1 columns of A[n], we get four linear endomorphisms 
p(s,t) £ End(w4), < s,t < 1, of the vector space A. We call these 
endomorphisms shift maps. Using a hopefully suggestive synthetic notation, 
an element A £ A can thus be written as 

A - A[ ° l \p(l,0)A p(lA)A 
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with A[0] G C and p(0, 0)A, p(0, p(l, 0)A, p(l, I) A G A 

Definition A subspace V C A is recursively closed if /9(s, i)V C V for all 

s, i. 

The recursive closure S of a subset S G .A is the smallest recursively 
closed subspace of A which contains S. We denote by A the recursive closure 
of the subset {A} reduced to a single element A £ A. The complexity of 
A G A is the dimension dim(A) G N U {00} of the recursive closure Ac A 

An element A £ A is a recurrence matrix if its recursive closure A is of 
finite dimension. We denote by 1Z C A the subset of all recurrence matrices. 

Writing p{X S)t )A or simply X Syt A for p(s,t)A, < s,i < 1, the shift 
maps p(s,t) G End(A) induce a linear representation (still denoted) p : 
{X ,o, Ao,i,Xi i0 ,^i,i}* — > End (.A), recursively defined by 

(^si,ti^s 2 ,t2 - - - X s u ti) A = ( X si,tiX S2jt2 ■ ■ ■ ^s;-iA-i)(p( s Zi U)A) , 

of the free non-commutative monoid {-Xo,o, -X"o,ij X^o, X^i}*, called the s/ii/i 
monoid, in four generators A^o, -X"o,ij -^1,0 > -^1,1 representing shift maps. 
Subrepresentations of p correspond to recursively closed subspaces V of A 
spanned by (unions of) orbits under {A"o ; o, -X"o,i> A^o, 

The linear action of the monoid {Xo ; o, Ao,i, -X"i,o> ^1,1}* on ^ suggests 
to consider the bijective map which associates an element A G A with the 
non-commutative formal power series 



X£{Xo,oA,li^i,Oi^i,i}* 

in four free non-commutative variables Ao,o, Ao,i, -X"i,0; This bijection 
restricts to a bijection between the vector space 1Z of recurrence matrices 
and rational elements in C((Xo ; o, Ao ; i, Xi ; o, A^i)). 

The algebraic structure of TZ C A is described by the following result. 

Proposition 3.1. (i) We have dim(XA) = dim(A) for all A G C* and for 
all A G A. 

(ii) We have dim(A + B) < dim(A) + dim{B) for all A,B G A. 
(Hi) We have dim(AB) < dim(A) dim{B) for all A,BeA. 

Remark 3.2. The inequalities of assertion (ii) and (Hi) can of course be 
strict: Consider two elements A, B G A defined by A[n] = Id[n], B[n] = 

1- ^2^ Id [re] where Id[n] denotes the identity matrix of size 2 n x 2 n . The 
elements A,B have common recursive closure A = B = CA + CB of dimen- 
sion 2. Their sum A + B = Id G 1Z is the identity element having complexity 
1 and their product AB = has complexity 0. 

Corollary 3.3. The set 1Z of recurrence matrices is a subalgebra of A. 
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Proof of Proposition 13.11 (i) and (ii) are obvious. 

Denoting (slightly abusively) by A B = {J2 x i Y i \ x i G ~A,Yi £ B} 
the vector space spanned by all products XY, X G A, Y £ B, we have 
AB eAB. 

For 

XY-(XMYM) ( ^°'°)( Xy ) P(0>1)(^) \ f in 

with leijFGB, the computation 

p(0,0)(XF) = (p(0 ) 0)X)(p(0,0)r) + (p(0,l)X)(p(l,0)F) 

P (o, i)(xy) = (p(o, o)x)( P (o, i)y) + ( P (o, i)x)(p(i, i)y ) 
/9 (i,o)(xy) = (p(i,o)x)(p(o,o)y) + ( P (i,i)x)(p(i,o)y) 

p(l, 1)(XF) = (p(l, 0)X)(p(0, l)y ) + (p(l, l)X)(p(l, l)Y) 

shows that A B is recursively closed of dimension < dim(A) dim(B) . As- 
sertion (iii) follows now from the obvious inclusion AB C A B. □ 

Remark 3.4. Certain properties of binomial coefficients (mod p d ) are easy 
to study using "recurrence matrices" given by sequences of matrices of size 
pi x p> , j = 0, 1, . . . with entries in the associative ring Z/p d Z. 

3.2 Recursive presentations 

An element A E A is completely determined by the action of the shift maps 
p(s,t) on its recursive closure A (spanned by {A^o, Xo,i, X^q, Xi t i}*A), 
together with the restriction to A of the augmentation map ttq G A* defined 
by projecting an element X = (X[0],X[1], . . . ) € A onto its initial value 
7ro(X[0],X[l],X[2],...)=X[0]eC. 

A recurrence matrix A can thus be given by a finite amount of data: 
An (expression describing the) element A of the finite-dimensional vec- 
tor space A, the restriction (still denoted) ttq £ A of the augmentation 
map expressing the initial values of elements in A, and a 2 x 2 matrix 

P = ( P ^?l I £ M 2 x2(A®A*) of tensors encoding the shift 

V Pi 1 * ) KM) / 
maps. The coefficients of the matrix A[n] S M2" X 2 n are then obtained by 

"contractions" of ttq p n A. 

This leads to the notion of recursive presentations. A recursive presen- 
tation for A E 1Z is given by the choice of a basis A\ = A, . . . , A a of A (or 
more generally of a finite set A±, . . . , A a spanning a recursively closed vector 
space containing A) and by recursive identities 

A = 7r (A \ ( PM A 3=T, a k=iP(0,0)k, j A k p(0,l)A j = ^2 a k=1 p(0,l) k , j A k 
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encoding the initial values iro(Aj) and the values p(s, t){Aj) G A of the shift 
maps. A recursive presentation defines the elements A% = A, . . . , A a span- 
ning (a recursively closed subspace containing) A recursively by expressing 
the four "blocks" of Aj [n + 1] as linear combinations of A\ [n] , . . . , A a [n] . 

Remark 3.5. We use the convention that G 1Z admits the empty presen- 
tation and is "the" element of an empty basis. We speak thus of "the" 
basis Ai, . . . of representing = A±. 

3.3 Saturation level 

We denote by tti(A) = A[l] the projection of a matrix sequence A = (A[0], A[l], A[2], 
A onto its square matrix A[l] of size 2 l x 2 l . Similarly, 

k<i(A) = (7ro(A),7ri(A), . . .,MA)) = (A[0],A[1],. . . , A[l}) g ® l j=0 M 2j x2i (C) 

denotes the projection of the sequence A onto its first / + 1 matrices. 

The saturation level of a finite dimensional subspace V C A is the small- 
est integer N G N such that /C<at(V) = IC<n+i(V) where /C<;(V) C V is the 
kernel of the projection 7r</ : V — > 7r<;(A) = ©j- = qA^2jx2j- 

Proposition 3.6. We /iaue /C<at(V) = {0} for the saturation level N of a 
finite- dimensional subspace V C A which is recursively closed. 
In particular, it<n : V — > ®^Lo-^2-?x2J defines an injection. 

Proof The obvious inclusions p(s, t)/C</ + i(V) C /C</(V) imply that /C<at(V) = 
/C<at_|_i(V) C V is recursively closed. Since the restriction to IC<n+i(V) C 
/C<o of the augmentation map -kq : A — ► C is trivial, we have (XA')[0] = 
for all X G {Xq^, Xq^, X\f h Xi t \}* and for all K G /C<jv(V). This shows 
JC< N (V) = {0}. ' □ 

Proposition 13.61 enables us to extract a basis from a finite set S spanning 
a recursively closed vector space V C 1Z. Similarly, Proposition 13.61 allows 
the construction of a basis of the subspace A C V for an element A G V of 
a finite-dimensional recursively closed vector space V C 1Z. 

These operations are the necessary ingredients for effective computations 
in the algebra 1Z. Effectivity means that there exists an algorithm involving 
only a finite number of elementary operations in the groundfield C and a 
finite amount of data which computes the result of an algebraic expression 
(given by a non-commutative polynomial) involving (recursive presentations 
of) a finite number of elements in 1Z. 

The necessary elementary algorithms can be briefly described as follows: 

3.3.1 Multiplication of A G 1Z by a non-zero scalar A G C* 

A presentation of XA is obtained from a presentation of A by multiplying 
the initial values Tr(Aj) = Aj[0] G C with A (and by keeping the same shift 
maps). 
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3.3.2 Addition of two elements A, B G 1Z 

Add a first element A\ + A2 having the obvious initial value tto(A + B) = 

Ax[0] + -Bi[0] to the list of not necessarily linearly independent elements 

Ax, ... , A a , Bi, . . . , Bb spanning A + B. The elements p(s, t)(A\+A2), p(s, t)Aj, p(s, t)Bj 

are given by 

a b 

p(s, i)(Ai +B 1 ) = Y, p(s, t)l % A k + p( s > t^iBk 
fc=i fc=i 

and 

a 6 

p(a,t)A = J2p( s > t )k,i A k, p{s,t)Bj = Y,p( s ^)kj B k 

k=l k=l 

for all < s, t < 1, 1 < i < a, 1 < j < 6 where p(s, t)" 4 and i) B are the 
obvious shift maps with respect to bases A = Ai, . . . , A a and B = B\, . . . , Bf, 
of A and B. Working with the finite-dimensional recursively closed subspace 
C(Ai + B\) + Yli=i CAj + X]j=i ^--^i c ^ one can now gi ve a presentation 
of Ai + B\ by computing a basis of A\ + B\, followed by the computation 
of the coefficients (with respect to this basis) of the shift maps. 

Remark 3.7. The algorithms VJ.'J. l\ and VJ.'J^ can be used to compare two 
elements A, B £ TZ by computing a presentation of A — B. 

3.3.3 Multiplication of two elements A, B G TZ 

Consider the ab elements Cj,- = AiBj, 1 < i < a, 1 < j < b with initial 
values Cjj[0] = no(AiBj) = Ai[0]Bj[0]. Shift maps are given by 

(*) P (s,t)c hJ = ZUZli (pMiA^j + pi^iM 1 ^ c ^ 

using the notations of 13.3,21 One constructs now a recursive presentation of 
C\ t i = AB by proceeding as above using the recursively closed vector space 

E"=i Ej=i 3 

Remark 3.8. The formulae (*) occuring in \3.3.3\ define an associative prod- 
uct on the set C of all equivalence-classes of finite- dimensional linear rep- 
resentations of the monoid {Xo^, ^o.i> ^i,o> ^l.i}*- Considering also di- 
rect sums of linear representations turns C into a semi-ring. The semi- 
ring C has a homomorphism ip into the semi-ring (with addition V + W = 
{X + Y I X G V, Y G W} and product VW = {£\ X { Y { \ X; G V, Y; G W}) 
formed by all finite- dimensional recursively closed subspaces oflZ. The ele- 
ments of ip(L) have two equivalent descriptions: They can be identified with 
the subset £ C C given by all equivalence classes of finite- dimensional linear 
representations of{Xo i o, Xq,!, -X"i,0) -^1,1}* containing no equivalence class of 
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a subrepresentation with multiplicity > 1. The second descriptions involves 
birecursively closed vector spaces which are defined as follows: A vector space 
V C 1Z is birecursively closed if it is recursively closed and if an arbitrary 
generic modification of the initial values in a recursive presentation of an 
element in V yields again a presentation of an element in V. 

The homomorphism of semi-rings tp associates to a finite- dimensional 
linear representation pj of {Xq^, Xo,i, Xip, Xi t i} the unique maximal bire- 
cursively closed subspace of 1Z whoses shift maps involve only equivalence 
classes of subrepresentations in pf. 

Finite- dimensional birecursively closed subspaces of 1Z are stable under 
direct sums and products and their semi-ring is the quotient semi-ring <p(C) 
ofC. 

3.4 The LU decomposition of a convergent non-singular ele- 
ment in 1Z 

An element P £ A such that P = p(0, 0)P is called convergent. It is given 
by considering the sequence 



P 



0,0) 



-Po.o -Fb.i 
Pi,o Pi,i 



( Pq,g Pq,i Pq,2 -Po,3 \ 

P\,Q P\,l Pl,2 Pl,3 

P2,0 P2,l P2,2 P2,3 

\ -^3,0 ^3,1 Ps,2 Ps,3 J 



of all square submatrices formed by the first 2 n rows and columns of an 
infinite "limit" matrix 

fo,o -fb.i Pq,2 ■ ■ 

Pl,0 Pl,l Pl,2 ■ ■ 



Henceforth we denote generally a convergent element in A and the associated 
infinite matrix by the same letter. This should not lead to confusions except 
in cases where both interpretations are correct. 

We call an infinite matrix P non-singular if the k X k square matrix P(k) 
formed by its first k rows and columns has non-zero determinant for all /c > 1. 
Such a non-singular matrix P has an LU— decomposition: It can be written 
as P = LU with L lower triangular unipotent (l's on the diagonal) and U 
upper triangular non-singular. The identity P = LU implies the equality 
det(P(fc)) = det(U(k)) for all k > 1 and gives rise to an LU— decomposition 
in A by considering as above for n = 0, 1, 2, . . . the submatrices formed by 
the first 2 n rows and columns of of P, L and U. If P is symmetric we have 
moreover U = DL l where D is diagonal non-singular and L* is obtained by 
transposing the matrix L. 

All proofs of the results presented in Section[2]boil down to LU— decompositions 
with P = P l = LU, L,D,U = DL l e 1Z. 
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Remark 3.9. Call an element A £ A non-singular if it involves only non- 
singular matrices A[0], A[l], .... Such an element has an LU — decomposition 
(in the obvious sense) in A. 

Proving the non-existence of an LU — decomposition in TZ for a suitable 
given non-singular recurrence matrix A £ TZ is probably difficult. 

The related problem of constructing the (existing) recurrence matrices 
L,U £ 7Z from the knowledge (of a recursive presentation) of A = LU £ TZ 
has however an algorithmic answer: One proceeds as for the existence of 
an inverse element by guessing recursive presentations for L and U using 
LU— decompositions of finitely many matrices A[0],A[1], . . . ,A[N + 1]. In 
case of succes, the resulting hypothetical decomposition, if correct, can then 
be proven to hold. The necessary algorithm is obtained after minor modifica- 
tions from the algorithm for computing the inverse A' 1 £ TZ of an invertible 
element A £ TZ described in Section\4-l\ 



4 Invertible recurrence matrices 

The set of all recurrence matrices which are invertible in the algebra TZ forms 
the group of units in TZ. Determining the inclusion in the unit group of TZ 
of a recurrence matrix A is perhaps a difficult problem without algorithmic 
solution. Indeed, we have the following result. 

Proposition 4.1. (i) For every natural integer n there exist invertible re- 
currence matrices A,B = A -1 £ TZ such that dim(A) = 2 and dim(I?) > n. 

(ii) There exist elements in TZ which are invertible in the algebra A but 
not in the subalgebra TZ of recurrence matrices. 

Remark 4.2. The assumption dim(A) =2 is optimal: Invertible recurrence 
matrices of complexity 1 form a subgroup (isomorphic to C* x GL2(C)) in 
TZ. 

Proof of Proposition 3H] For uj £ C*, consider the convergent element 



/ 


1 




\ 




— u 


1 









—UJ 


1 


V 








-UJ 1 J 



, • • • £ A 



consisting of lower triangular unipotent matrices with constant subdiagonal 
—uj. It defines a recurrence matrix A = A± of complexity 2 recursively 
presented by 

, / M \ . / A 2 \ 
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Since A is given by a sequence of unipotent lower triangular matrices, it 
is invertible in the algebra A with inverse the convergent element 



/ i 



UJ 
2 



V UJ 3 



UJ 
UJ 2 



G A 



UJ I J 



whose limit is the infinite unipotent lower triangular Toeplitz matrix with 
constant subdiagonals associated to the geometric progression 1, uj, uj 2 , . . . . 
For k a strictly positive integer, we consider the element 



, ,2°k , 2 1 k 
UJ , UJ 



1 



UJ 
2 



2 2 fc 



/ 


1 


UJ 


uj 2 


UJ 3 




UJ 


UJ 2 


UJ 3 






uj 2 


UJ 3 


uj 4 


UJ 5 


V 


UJ 3 


UJ 4 


UJ 5 


UJ G 



\ 



G A 



u 2 n k+i+j^ < i, j < 2™. A straightforward computa- 
tion shows p(s,t)Sk = S s+ t+2k for all s,t. Since we have 



given by (S k [n]) id 



B 



B 
Si £ 



(using the notation of recursive presentations) , the algebraic closure B of B 
is spanned by the set 

{X 0fi ,X 0jl ,X lfi ,X 1A }*B = BU {S k | k> 1} . 

Since a;^ = 1 implies Sfc+^r = S&, this set contains at most 1 + elements 
if uj is a root of 1 of finite order N. 

In order to prove assertion (i), we consider the case where w is a root 
of 1 having odd order N > 2 n and we denote by a > log 2 (N + 1) > n the 
order of 2 in the multiplicative group (Z/iVZ)*. The action of the monoid 
p(0, 0) N on Sk corresponds then to the action of the Galois map uj i — ► uj 2 
on the upper left coefficients 

, 2°fc , 2!fe , ,2 2 fc , ,2 3 fc 

CJ , UJ ,UJ ,UJ , • • • 

of Sk- Elementary number theory (for instance reduction modulo 2 by choos- 
ing uj among the primitive N— th roots of 1 in the field extension ¥ 2 a of de- 
gree a over F2) shows that the set S 2 n spans a subspace of dimension a > n 
in B. This implies assertion (i). 

Consider now A as above with wGC\Q transcendental. This allows to 
consider a; as a variable and we have dim(S^) > dim(B^) for the complexity 
of the inverse element B w = A^ 1 where £ € Q is any algebraic specialisation 
of uj. Assertion (ii) follows now from assertion (i). □ 
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Remark 4.3. An example of B = A~ x 6 A \ 71 with A G 1Z invert- 
ible only in A is also given by A, B as above with uj = 2 (the argu- 
ment below works in fact for any uj of norm \uj\ > 1). Indeed, other- 
wise, up to a constant, the initial values 2, 2 2 , 2 4 , 2 8 , . . . of the sequence 
p(l,0)B, p(0,0)p(l,0)B, p(0,0) 2 p(l,0)B, . . . should be bounded above by a 
geometric progression l,p,p 2 ,... for any positive p exceeding the spectral 
radius (absolute value of the largest eigenvalue) of p(0, 0) G End(i3). 

The simplest element of 1Z with an inverse in A\TZ is perhaps given by 
the central diagonal recurrence matrix A of complexity 2 defined by A[n] = 
[n + l)Id[n] where Id[n] denotes the identity matrix of size 2™ x 2 n . We give 
two proofs that the inverse element A" 1 , given by A _1 [n] = -^^Id[n], has 
infinite complexity. 

A first proof follows from the observation that the sequences of diagonal 
coefficients n+ \ +1 = (p(0, 0) fc ^4 _1 )[n] form the rows of the Hilbert matrix of 
infinite rank with coefficients Hij = , < 

In order to give a second proof, we observe that all coefficients of A^ 1 
are rational numbers and every prime number appears as a denominator in 
a suitable coefficient of A -1 . This is impossible for an element B 6 1Z with 
rational coefficients. Indeed, such an element B has a recursive presentation 
with data (consisting of initial values and coefficients of shift maps with 
respect to a basis C {Xq$, -Xo,i, X\^, Xi^}*B of B) given by a finite set V 
of rational numbers. Since all coefficients of B are evaluations of integral 
polynomials on V, all denominators of coefficients in B involve only prime 
numbers occuring in the denominators of the finite set V. 

Remark 4.4. The quotient group of invertible elements modulo C*Id can 
be turned into a metric group by considering the positive real function 

A i — HI A 11= max(log(dim(]4 + C/d)),log(dim(]F T + C/cO)) 

on the group T of units in TZ. It satisfies \\ A \\> with \\ A ||= only 
for A G C*Id, || AB \\<\\ A || + || B || and \\ A \\ = \\ A' 1 \\ and defines 
thus a left-invariant distance on the quotient group T/C*Id by considering 
d(A, B) =|| A- X B ||. 

The corresponding group over a finite field has finitely many elements in 
balls of finite radii and it would be interesting to understand the generating 
function 



where the sum is over all elements of the unit group r(F pe ) of the algebra 




t 




A£T(F p e) 



TZ(¥ p e) defined in the obvious way over the finite field F p e. 
The related generating function 



£ t dim(A) £ N[[t}] 



A€TZ(¥ p e) 
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counting all elements of given complexity in the algebra IZ(¥ p e) is probably 
fairly easy to compute. It has convergency radius and is thus transcenden- 
tal. 

4.1 Algorithm for computing ^4 1 G 7Z for A G 1Z invertible in 

K 

The following algorithm computes the inverse of an element A £ 1Z if it has 
an inverse in 1Z and fails (does never stop and uses more and more memory) 
for an element A as in assertion (ii) of Proposition 14.11 Non-invertibility in 
A will eventually be detected (assuming exact arithmetics over the ground 
field) by exhibiting an integer k for which A [k] is singular. 
Input: A (presentation of a) recurrence matrix A G 1Z. 

Set N = 1 and M = 1. 
Loop Compute, if possible, the matrices B[0] = Aft]- 1 , B[l) = A[l] _1 , . . . ,B[N+ 
M + l] = A[N + M + 

If a matrix A[k] with k < N + M + 1 is not invertible, print "The matrix 
A[k] is singular and A has thus no inverse in .A" and stop. 

We denote by B G A the sequence B[0], . . . , B[N + M + 1] completed 
by arbitrary matrices B[k] of size 2 k x 2 k (for example by zero matrices) if 
k > N + M + 2. 

For n, m G N such that n + m < N + M + 1, we denote by V(n, m) C 
ffij = o-^2ix2J the vector-space spanned by all 1 + 4 + • • • + 4 m elements 
(vr< n (XB)) Xe ^< m where X- m denotes the set of all = 1 + 4 + 4 2 + 

..._)_ 4 m words of length < m in the alphabet ^ = {-Xo,o, -^0,1 > -^l,o> ^1,1}- 

If dim(V(iV, M)) < dim(V(iV + 1, M)), then increase N by 1 and iterate 
the Loop. 

If dim(V(A r , M)) < dim(V(iV, M + 1)), then increase M by 1 and iterate 
the Loop. 

Setting d = dim(V(iV, M)), the identity dim(V(iV,M)) = dim(V(iV + 
1, M)) shows that the finite sequence B[0], . . . , B[N + M + 1] can be com- 
pleted to a uniquely defined element B G 1Z which is of complexity d and of 
saturation level < N . Use the natural isomorphisms between V(N, M), V(iV+ 
1, M) and V{N, M+l) and the inclusions p(s, t)V(N, M+l) C V(N, M), < 
s,t < 1 (where p(s, t) acts in the obvious way) for writing down a recursive 
presentation of B. 

Check if AB = 1 using the algorithms of Section 13.31 

If yes, print the presentation found for B and stop. 

Otherwise, increase N by 1 and iterate the Loop. 
End of Loop 

Remark 4.5. The above algorithm can perhaps be improved. In particular, 
it is probably not necessary to consider all (4 M+1 - 1)/3 words of X^ M C X* 
during the computation of the dimension ofV(N + 1,M). 
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5 Modulo 2 and 2— valuations 



Proof of Theorem 12.11 The infinite symmetric Pascal matrix P with co- 
efficients (^(*^) (mod 2)^ G {0, 1} for < i,j defines a convergent element 

(still denoted) P £ A. It follows from Lucas's formula (see Section T2.ip that 
P is a recurrence matrix of complexity 1 recursively presented by 

f P P 

(zero-entries are omitted). The recurrence matrix P G 7Z has an LU decom- 
position in TZ given by the equality P = LDL 1 with L, D G TZ of complexity 
1 defined by the recursive presentations 

L = 1,( L L L ) wLD = l{ D _ D 

where L is lower triangular unipotent and D is diagonal. An easy analysis 
of the coefficients of the diagonal recurrence matrix D ends the proof. □ 

Remark 5.1. The convergent lower triangular recurrence matrix L and the 
convergent diagonal recurrence matrix D correspond to the infinite limit- 
matrices (still denoted) L,D with coefficients given by L{j = (mod 2)^ 
{0,1} (for < i,j) and D n>n = (_l)^o+^i+^2+... _ (_x)^(n) where n _ 
^2j>o l/ j^' > ^ is a binary integer. 

Remark 5.2. Recurrence matrices of complexity 1 are, up to scalars, of the 
form 1, M, M (8) M, M ® M ® M, . . . where M is a complex 2x2 matrix. 

It follows that the matrices L and D ( and thus also P = LDL 1 ) involved 
in the proof of Theorem \2.1\ are invertible in TZ. The recurrence matrix 
D is its own inverse. The inverse L^ 1 of L is recursively presented by 
L- 1 



1 - L] I. I. 1 

Remark 5.3. The spectrum of recurrence matrices of complexity 1 is easy 
to compute: Given a square matrix M of size d x d with characteristic 
polynomial YYj=i (^ — ^j) we have 

n (*-**!••• 

(h,...,j n )e{i,...,d} n 

for the characteristic polynomial of the iterated tensor power M® n . 
Proof of Theorem 12.31 The infinite matrix V with coefficients 

y st = ^2(( S t')) = jds(s)+ds(t)-ds( S +t) 
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gives rise to a convergent element V £ A. A bit of work using Kummer's 
formulae (see Section [2. 2p shows that V = V% is recurrence matrix given by 
the recursive presentation 

Vl - l \v 2 iVt 

v = 1 ( Vi -iVi + (1 + i)V 2 

2 ' I -^l + (1 + ^2 -Vi 



We have V = LDL 1 G TZ with L = L\ £ TZ recursively presented by 
Li = 1, 



L 2 = 0, 
L 3 = l, 
L 4 = l, 



Li 

L3 L4 

-iL 2 
— L\ + L 3 — iL 2 — iL± 

L\ L 2 
-iLi + (1 + i)L 3 L 2 + (l + i)L 4 

(1 - i)Li + fL 3 L 4 



and with diagonal D = D\ G TZ recursively presented by 



L>i = 1 

D 2 = -l + i, 
D 3 = -l + i, 



^2 

£3 

2D! - D 2 + 2D 3 

£3 

-D 2 



An analysis (left to the reader) of the diagonal entries of D\ ends the proof. 
□ 

Remark 5.4. The recurrence matrices L,D and V = LDL 1 are invertible 
in TZ, see J3]/. 

6 Beeblebrox reduction 

This section is devoted to proofs and complements involving the Beeblerox 
reduction Xb((T\) of binomial coefficients. 
Proof of Theorem 12.51 We have 

2n\ (2rc) • • • (2n - 2k + 1) (n\ (2n - l)(2n - 3) • • • (2n - 2k + 1) 



K 2kJ {2k) ---I \kj (2fc-l)(2fc-3)---l 

where both the numerator and the denominator of the fraction 

_ (2n - l)(2n - 3) • • • (2n - 2k + 1) 
~~ (2fc — l)(2fe — 3) ■ ■ • 1 
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contain k terms. If k is even, we have F = 1 (mod 4) since the numerator 
and denominator of the fraction F contain both k/2 factors = 1 (mod 4) 
and k/2 factors = — 1 (mod 4). If A; and n are both odd, the numerator and 
denominator of F contain both (k + 1) /2 factors = 1 (mod 4) and (k — l)/2 
factors = — 1 (mod 4) and we have again F = 1 (mod 4). If k is odd and n 
is even, then both binomial coefficients Qj) an d ()!) are even and we have 
thus Xb((^1)) = Xs((fc)) = 0- This proves the first equality. 

The binomial coefficient { 2 k+i) is obviously even and this implies the 
second equality. 

In the next case we have 

(2n + 1\ _ /n\ (2n + l)(2n - 1) • • • (2n - 2k + 3) 
V 2k )~\k) (2fc-l)(2fc-3)---(l) 

and the last fraction equals 1 (mod 4) if k is even. 

For k odd and n even, we have xM^fe 1 )) = Xb{(X)) = since f 2 ^ ) = 
©=0 (mod 2). 

For n, A; both odd, the correction (— l) fc = —1 equals the fraction modulo 
4. This ends the proof of the third equality. 
In the case of the last equality, we have 

(2n + 1\ _ 2n + 1 /2n\ 
\2k + lj ~ 2fc + 1 V2fc y 

which is even if n = (mod 2) and A; = 1 (mod 2). If n = k (mod 2) then 
§f±± = 1 (mod 4). For n odd and k even we have = -1 = (-l)"^ 1 ) 
(mod 4). The first equality and these observations complete the proof. □ 

6.1 Proof of Theorem 12.61 

Proof As in Section 13.41 we consider the element (still denoted) Z £ A 
associated to the infinite matrix Z with coefficients Z s ^ = Xs(( s ^*))> < 
s, t. 

We have to show that Z is a recurrence matrix and we have to find 
a recursive presentation for Z. This can be done in the following way: 
We consider left shift-maps A(0, 0), A(0, 1), A(l, 0), A(l, 1) which associate to 
A = (A[0],A[1],...) G A the element X(s,t)A G A where (X(s,t)A)[n] is 
the submatrix of A[n + 1] corresponding to row- indices = s (mod 2) and 
column-indices = t (mod 2). A subspace V C A is left-recursively closed 
if it is invariant under all four left shift-maps and the left-recursive closure 
A C A of A £ A is the smallest left-recursively closed subspace containing 
A. Theorem 12.51 implies that Z is finite-dimensional. One shows that 
dim(A ) = dim(A) for all A G A (see for example [J| for the details). 
This proves that Z is a recurrence matrix. A little bit of work based on 
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properties of the saturation index shows now that Z = Z\ is given by the 
recursive presentation 



We have the identity Z = LDL} with L = L\ £ TZ lower triangular 
unipotent given by the recursive presentation 



An easy inspection of the diagonal entries of D\ completes the proof. □ 

Remark 6.1. The birecursively closed subspaces of A appearing in Remark 
\3.8\ are the recursively closed subspaces of A which are also left-recursively 
closed, ie. invariant under all four left shift-maps mentionned above. 

6.1.1 The group T L 

All four recurrence matrices Li, L3, L4 involved in our recursive presen- 
tation of L = L\ are invertible in 1Z. Their inverses are given by 





The diagonal matrix D = D\ G 1Z has recursive presentation 





-1 



Mi, L 2 



-1 




1 



-M 2 , LJ 1 = M 4 



with Mi, M2, M3, M4 recursively presented by 




We have the curious inclusions 



p{0,0)L£QL 1 
p(l,0)LeQL 3 



p(0,l)LeQL 2 , 
p(l, 1)L G Q-L4 
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for L G L\ = L2 = 
for their inverses 

Mi 



L3 = L4 



:1 CLj. The analogous property holds also 



Mo = ~2Ln 



M, 



-2L.7 1 , M 4 = L7 1 . 



This suggests that it would perhaps be interesting to understand the 
group Tl = (a, b, c,d) C K generated by the normalised recurrence matrices 



b=-L 2 



1, 



1 

1 1 



1 

1 

3 


1 



£3 = 1, 



d = L4 = 1, 



V 

/ 1 

1 
1 
1 



\ 



1 / 

-I 4 \ 

I 2 ' 

1 
\ 1 



/ 1 
1 

v! 



4 \ 

2 





1 / 



corresponding to the invertible "projective" elements Q*Lj. In particular, 
it would be interesting to understand if Tl is a linear group. This would 
certainly be implied by the existence (which I ignore) of a natural integer 
N for which the projection Tl — ► vr<Ar(r^) is one-to-one. 

L. Bartholdi communicated to me the following list implying all relations 
of length < 12 in T L : 



ca 

~ca~ 
-1, 



-1. 



bd~ 
ca' 

cad~ 1 a~ 1 d a~ 1 b~ 1 , 
cad~ 1 c~ 1 bda~ 1 b~ 1 , 
cda~ 2 dad~ l b~ l , 
ca 2 dT 1 aT 2 dada~ 2 b~ 
d 2 ad~ 



l b 1 cada 3 , 
cdadT 7 'a^ 1 dada~ 2 b~ l . 

He observed that they are all of the form u n v~ l w n x~ l , where u n , v n , w n , x n 
are positive words of length n with respect to the generators {a, b, c, d}. 

More generally, it should also be interesting to understand the subalgebra 
C C 1Z generated by the recurrence matrices a 1 * 11 , 6 ±1 , c ±1 , d ±1 G 1Z. 
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The algebra C is of course a quotient of the group algebra C[Tl] and it 
would be interesting to describe the kernel of the associated homomorphism. 
The algebra £ is a quotient of the free non-commutative algebra 

M = C(A, A~\B, B~\C, C~\D, D- 1 ) 

in eight free non-commutative variables A ±l , B ±l , C ±1 , D^ 1 . The corre- 
sponding natural homomorphism it : M — > C (given by Z ±x i — > z ±x G 1Z 
for (Z,z) £ {(A,a),(B,b)(C,c),(D,d)}) factorises through the group alge- 
bra of the abstract group Tl- The kernel X = ker(-7r) C N contains thus all 
relation of Tl and in particular the trivial relations 

AA~ l - 1, BB- 1 - 1, CC~ l - 1, DD- 1 - 1 . 

In order to gain some information on X, we can consider the morphism 
of algebras p\ : N — ► A^2x2(AT) given by 

A \ j ( A' 1 



! c I) )■ 1 ' V r> ■ i) 1 

B ~ ( -i/„ 3 2 £ ^ , B-i ^ - ( 2C " 



C/3 D J ' I D-73 



, ^ 2B\ 1 / -A- 1 2C- 1 

s u 1 C D J ' G ' V. B 1 D ; 



C/3 D J' \ -5-73 Z)- 1 J ' 

This morphism factors through it and induces a homomorphism of alge- 
bras ~pi : C — ► M2x2(£) which removes simply the first matrix X[0] from 
an element X[0], X[l], • • • £ C. This is due to the definition of fii which 
corresponds to the maps 1Z — ► M2x2(7Z) given by 



X 



p(0,0)X p(0,l)X 
p(l,0)X p(l,l)X 



for X £ {a ±1 , 6 ±_1 ,c ±1 ,d ±1 }. We have thus in particular pi(2) C M 2X 2(X). 
Application of p\ to some known element in X can sometimes be used for 
the discovery of new elements in X: The computation 

p x {AA- 1 ) = ( ^ D ^ ^ = ( CA A! £,5-1 £,£)-! ) 

shows that application of p\ to the trivial relation AA~ l — 1 € X implies the 
already known inclusions AA~ X — \,DD~ l — 1 G X and produces the non- 
trivial relation CA~ X — DB~ l £ X. Since elements of X involving only two 
monomials induce relations on the group Tl, we get the relation bd~ 1 ca~ 1 
in Tl which is the first relation in Bartholdi's list. 
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"Iterating" the map \i\ produces homomorphisms \i n : M — ► M 2 n x2 n(N) 
with similar properties. In particular, we have /J, n (I) C M 2 n x2 n(I). Index- 
ing the coefficients of matrices inM2" X 2 n (AA) by elements X £ {^o,o> -Xo,i> -^i,o> Xii} n , 
we get linear maps /ix : — > J\f by considering the coefficient correspond- 
ing to X in n n (Af) . The reader should be warned that the map X i — ► /ix G 
End{M) is not a morphism of monoids from {Xo 5 o, ^o,i> ^i,o> ^1,1}* m to 
End(M). The monoid generated by all maps fix, X G {^o,o> -^o,i> -^1,0 > -^1,1}*) 
preserves however the ideal T. 

Denoting by JV(Z) the subring of noncommutative polynomials with in- 
tegral coefficients, relations of the form r\ = r 2 in are in bijection with 
pairs of roots in the infinite-dimensional Euclidean lattice (with respect to 
the orthonormal basis given by monomials) In JV(Z). 

The peculiar form of all relations in Bartholdi's list is partially explained 
by the formulae for fi±. They imply that the maps /ix preserve sign struc- 
tures: The vector space spanned by the orbit under the monoid generated 
by the maps /ix of a relation of the form UnV^ 1 = XnW^ 1 contains only 
relations of the same form. 

It would be interesting to know if the ideal 1 is finitely generated as 
an {/ix}xe{x ,x i,Xi ,x 1 1}* —module: Otherwise stated, does I contain a 
finite subset Q such that X is the smallest bilateral ideal which contains Q 
and which is preserved by all maps /ix? 

Remark 6.2. The techniques used in this Section can of course be applied to 
other subsets in 1Z. One needs a (preferably finitely generated) algebra S (eg. 
the algebra generated by suitable elements of a subgroup inlZ) such that the 
recursive closure of every element in S is spanned by elements of S. The 
choice of a generating set Q allows to consider the free non- commutative 
algebra N on Q which gives rise to the natural surjective homomorphism 
7r : J\f — ► S . Choosing lifts of the shift maps, one constructs homomor- 
phisms /xpj : N — ► M 2 n x2 n(J\[) giving rise to the linear maps /ix £ End{M) 
preserving the bilateral ideal T = ker(7r). 

The maps /ix can be choosen in order to preserve the grading of Af if 
the generating set Q spans a recursively closed subspace. 



6.1.2 The inverses of Di,D 2 and the group Tz 

The inverses of the diagonal recurrence matrices D\ , D2 (defined by the 
decomposition Z = LD\L l and by D 2 = p(l,l)Di) are E\ = -Df 1 , E 2 = 
D^ 1 recursively presented by 

E, = l,( El „ ),E 2 = -lJ & 



E 2 J ' \ 3E 2 

The inverses of the matrices Zi,Z 2 ,Z% are U\ = Z-f 1 , C/3 = Z^ 1 ,U 2 
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Z^ 1 recursively presented by 

U 2 



' 1 1 u 3 u 1 -u 2 -u 3 



u 2 = i, 
u 3 = i, 



u 2 
-u 3 -u 1 + u 2 + u 3 
-u 2 
u 3 -u l + u 2 + u 3 



Since Z\ = and Z 2 = Z|, the group Tz = (a,b,c) C TZ generated by 
a = Z\,b = Z 2 ,c = Z 3 has an involutive automorphism given by a i— ► 
a -1 , 5 i — > c , c i — ^ b . Two relations in Tz are (ab -1 ) 2 and ab = ca. 
Using the relation a 2 = cb following from the computation 



a: 



2 — aab l b = aba l b = caa 1 b = cb , 



every element of Tz can be expressed as an element of a e (b, c) with e £ {0, 1}. 

Since det(ir 2 (a)) = — 1 and det(7T2(b)) = det(7T2(c)) = 1, the subgroup 
generated by b, c is of index 2 in Tz- 



6.2 The triangular Beeblebrox matrix 

We define the lower triangular Beeblebrox matrix as the infinite lower tri- 
angular matrix with coefficients L s j = XB(( s t )),0 < s,t, given by the Bee- 
blebrox reduction of binomial coefficients. 

One of the main results of [7] states that any fixed row of L contains 
either no coefficients —1 or the same number (given by a power of 2) of 
coefficients 1 and —1. This can of course also be deduced from Theorem 
12.51 or by computing L J where J is the "recurrence vector" obtained by 
considering the sequence of column vectors 

(1), (1,1)', (1,1,1,1)', (1,1,1,1,1,1,1,1)',... . 

The triangular Beeblebrox matrix L defines a recurrence matrix (still 
denoted) L = L\ 6 TZ recursively presented by 

6.2.1 The recurrence matrices L^ 1 

The lower triangular recurrence matrices L\,L 2 ,L 3 defined above are invert- 
ible in TZ with inverse elements Mi = L^ 1 , M 2 = , M 3 = LJ 1 recursively 
presented by 

Mi = l,(^ „ ),M 2 = 1,( M n l „ ),M 3 = 1,( ™] ° ) 
1 ' \ M M 3 J ' 2 '\ -M -M 3 / V —M M 3 J 

where M = M 1 - M 2 - M 3 . 
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Proposition 6.3. The map 

( 1 2 \ / 1 1 \ / 1 \ 

L\ i — ► J 1 ,L 2 ^ 1 1 \ ,L 3 ^ 1 2 
\ 1 / V 1 / \ 1 / 

(where the three matrices correspond to the three affine maps (x,y) *— > (x + 
2, y), (x, y) i — ^ (y + 1, x + 1), (x, y) (x, y + 2) o/M 2 j defines a faithful linear 
representation of the group {L\, L2, L%) C 1Z generated by L\, L2, L3. 

Moreover, the group homomorphism Li 1 — ► ^(Lj) is faithful on (L\, L2, L3). 

Proof We check that L\ and L3 commute. They generate thus an abelian 
subs; roup T a which is easily seen to be free abelian of rank 2 by considering 
the 4x4 matrices ^(Li) and ^2(^3). Checking the relations 

L\ = L1L3, L2L3 = L1L2, L2L1 = L3L2 

shows that r a is of index 2 in (L\, L2, L3) and these relations define the 
affine group of Proposition 16.31 

Faithfulness of the homomorphism Lj 1 — ► 7T 2 (ivi) follows from the obser- 
vation that the subgroup generated by ^2(^1)^2(^3) is free abelian of rank 
2 and does not contain ^2(^2)- a 



7 On the Jacobi-Dirichlet character 

This section contains the most important data for proving Theorem 12 .71 We 
omit the somewhat lengthy details. 

Tedious work proving formulae analogous to Theorem 12.51 or general 
principles show that the matrix J = J\ with coefficients Xj(( s ^'))' ^ — s ^ 
is of complexity 9 and has a recursive presentation given by 

*--(3*)- 

■<> ■•( J ')- ■<'< ••(£ Jl )- 

•'• '•( % -'")■ ■■■{$ '"■)■ 

where X 1 = X[0]*,X[1]*, X{2}\ ... for X £ A. 

The matrix J has an J = LDL 1 decomposition in 1Z with L of complexity 
20 and D of complexity 4. 
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The lower triangular recurrence matrix L = L\ involved in the decom- 
position J = LDL 1 has the recursive presentation: 



L 3 = l,| L " ° I, U 

Lo, L\q 


/. l.( Ls ; Ll2 L T 15 i i i. 

Lo, L\q 
1/3 '( 1/3L 6 )' Ll0 = 1 ' 

Ll6 ^13 



r r 

L4 L$ 






), 

/ 


L\ L\\ 


V 


L2 L 3 




La L14 




Lq L 7 


)• 


Li 





1/3L 2 


L3 j 


L 8 





1/3L 9 


L10 


f L19 


4L 5 


V 






L U = 2, 7 iD r w , L 12 =4/3, 
' ^17 -^18 / 

2/3Li2 2Li3 \ / L12 + Li6 L20 

; !! J' 7 ' ;; /-.v Lis 

, .-, ; — 2/3Li2 + Li6 L13 \ r o/q 1 — 2/3Li 2Ln 

^15 = A I 1 /nr r , -^16 — — ^/•J, 



L17 


= 2/3, ( 


L18 


-•(" 


L19 


= 8/3, ( 


L20 





1/3L 17 Lis / ' ' ' V 2 /3L 2 2L 3 

-2/3L 4 -4L 5 + 2L M 
2/3L 6 2L 7 
-2/3Lg — 2/3Li2 2L15 
2/3Lg 2Lio 
2L19 8/3L5 


4/3L 12 - 2/3L 16 -2Li3 + 2L 20 
2/3Li7 2Lig 

The four matrices L5, Li 2 , L13, L19 span a somewhat trivial four-dimensional 
subalgebra consisting only of matrix sequences with zero coefficients except 
for the first row. 

Consideration of the images p(s, t)L yields the following decomposition 
of the recursively closed vector space L = 0| =1 CLj: 

p(0, 0)V = CLi © CL 4 © CL 8 © CL12 © CLi 6 © CL19 , 
p(0, 1)V = CL 5 © CLn © CL13 © CL M © CL15 © CL 20 , 
p(l, 0)V = CL 2 © CL 6 © CL 9 © CLi 7 , 
p(l, 1)V = CL 3 © CL 7 © CL10 © CLis . 

I ignore if the vector space L contains generators of interesting algebras 
or groups. 

The diagonal recurrence matrix D = D\ involved in J = LDL 1 is of 



27 



complexity 4 with recursive presentation 



D 2 ) ' ' 'V Da 

8 binomials 

Proof of Theorem 12.91 The result holds for b = or for a < b. An 
induction on a + b ends the proof. It splits into the four following subcases: 
If a, b (mod n): 

a), = " 6 (VL+(ri\ 

_ /[a/nh / 6/a-l (mod n)\ . /'a— 1 (mod n)\ \ 

~~ UV«J^ V V b (modn) ^"""Vb-l (mod n)/ w j 

(\a/n\\/a (mod n)\ 

~~ \[b/n\)\b (modnj/^ 

If a = (mod n),b ^ (mod n): 

® q = " 6 (VL+(r;i 

_ (a/n— 1\ / {,/ n-1 \ . / n— 1 \ \ 

~~ \[b/n\ ) \ ^ \b (mod n)^ + 1(6-1) (mod n)) u J 

_ /a/n— 1\ I n \ 

~ \ [b/nj ) \b (mod n)) u 

an< ^ (& (mod n)) = ^ since 6^0 (mod n) implies that it is divisible by the 
n— th cyclotomic polynomial. 
If a ^0,6 = (modn): 



/a\ b(a— 1\ i ^a— 1 



fa— 1"^ , / L a / n J \ (a—1 (mod n)^ 
+ \b/n-l) \ 



b /ui ^ \b/n-l) \ n-1 

i/r, 
i/r. 
i/r 
b/r. 



([a/n\^a-l (mod n)^ _j_ q 



/La/"-J\ /a (mod n)\ 
\ b/n )\ J 



If a = 6 = (mod n): 



/a/n— 1\ /n— 1\ , (a/n— 1\ (n—l\ 

\ b/n )\ Q ~ l ~ U/n-lAn-l/ t 

i/r. 
U/r, 



/a/n\ /0\ 



□ 



Remark 8.1. / thank C. Krattenthaler for pointing out that Theorem \2.9\ 
follows also from the classical q— binomial identity 



n—1 n / v 

n(w<)=EH)Vs>n . 

,--n J,-n W o 



i=o fe=o 
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Setting n = a and q = u for lo a primitive d—th root of 1 and equating the 
coefficients of t b on both sides implies the result easily for odd d. The case 
of d even requires also a sign analysis. 

The above identity follows by induction on n from the easy computation 



n- =0 (W*) = (i-^)ELoH) fc <7 



EJ2(-*)V 



Theorem 12.91 implies in particular that for u = e 2lnk / n ^ the complex 
numbers (f)^, a, b G N belong to the finite subset U < b < a<n R> (l ( ( ™°^ "j) 
of real half- lines in C. 

For n = 2, the matrix (£)_., , < a, b < 2 is given by 



1 
1 1 



This implies that G N for all a, b G N. 

For n = 4 we get (£) . G N U iN U (1 + i)N since we have 

/ 1 \ 
1 1 

1 1 + i 1 

V i % % i / 



for the matrix with coefficients < a, 5 < 4. 



8.1 Tensor products 

Replacing coefficients in the ground-field C by multiples of a fixed matrix 
X of size 2 x 2^, one can consider the tensor product 

A®X = (A[0] X, A[l] ® X, . . . ) 

of yl G .4 or A G 1Z with X. Such an element has an LU decomposition 
A® X = {V (8) Lx){U' <S> Ux) involving elements of the same form if and 
only if we have decompositions A = L'U' G 1Z and X = LxUx- 

Remark 8.2. More generally, one can consider the quotient algebras A/J-S 
and 1Z/J-S where J-S is the ideal of all matrix sequences in A which in- 
volve only finitely many non-zero matrices. Elements of the quotient alge- 
bra IZjTS can be represented as linear combinations of suitable elements 
in Uk^IZ <S> M 2 k x2 if and such representations are sometimes simpler than 
recursive presentations of preimages in 1Z. 
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8.2 The Beeblebrox reduction of ('^ t )_ 1 

We denote by Z' the infinite symmetric matrix with coefficients Z' s t = 
Xs(( S ^')-i) §i ven by the Beeblebrox reduction of q— binomials evaluated 
at q = — 1. 

Theorem and Section O imply that Z' = L'D'iL'f where 

Z' = Z®(^ l l 1 ) = L®( 1 1 A, D' = D® 

(the tensor product X <g> M denotes the matrix(-sequence) obtained by re- 
placing a scalar entry A of Z by the 2x2 matrix M) with Z, L, D as in the 
proof of Theorem 12.61 

In particular, using Theorem l2.6l and the LU decomposition of 

( XB(( s+t ) 1 ) ) , one can easily write down a formula for det(Z'(n)) G 

V v s /_i 7o<s,{<l 

±3 N with Z'{n) denoting the symmetric n x n submatrix consisting of the 
first n rows and columns of Z'. 

Remark 8.3. The case of the matrix (with coefficients in {0, 1}) obtained by 
reducing modulo 2 yields nothing new since = = ( S ^*) 

(mod 2). 

8.3 Reduction modulo 2 and Beeblebrox reduction of ( s +*) . 
Let M' be the symmetric matrix with coefficients 

^(( S ^) ) €{0,±l,±x,±y}, 0<s,t 

where 

iP(^) = } j(a)x if £ = ai G iN 

( 7 (a)y if £ = a(l + ») G (l+t)N 

where 7 : N — * {0, ±1} is either the reduction modulo 2 with values in 
{0, 1} or the Beeblebrox reduction \B- 
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As in section E2] we have M' = L'D'iL') 1 where 



M' = M 



V 



D' = D 



(1 1 1 1\ 

1 y x 
1 x 



\ 1 

/ 1 
1 

1 

V i 

/ 1 



V 



/ 



1 

1-x 

l-y 5 

y-i 



l 



-2x+y 



1 / 



x — 2x+y 

i-y 



-2x+y 



J 



M, L, D are given by the matrices P, L, D, respectively Z, L, D occuring in 
the proof of Theorem 12.11 respectively 12.61 if 7 is the reduction modulo 2, 
respectively the Beeblebrox reduction. 

It follows that the determinant det (M'(n)) of the finite matrix M'(n) 
consisting of the first n rows and columns of M' is of the form 



±3 N x 2N (y-l) 



{o,i} (x 2_ 2x + y){ o,i} 



with powers of 3 only involved if 7 is the Beeblebrox reduction. The factor 
(y — 1) appears if and only if n = 2 (mod 4) and the factor (x 2 — 2x + y) 
appears if and only if n = 3 (mod 4). Using Theorem 12.11 respectively 
Theorem 12.61 it is easy to write down a formula for det (M'(n)). 

Acknowledgements I thank L. Bartholdi, M. Brion, P. de la Harpe 
and C. Krattenthaler for helpful comments and discussions. 
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